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MODELING OF COMBUSION PROCESSES OF STICK PROPELLANTS 
VIA COMBINED EULERIAN-LAGRANGIAN APPROACH. 


K. K. Kuo, K. C. Hsieh, M. M. Athavale 
Department of Mechanical Engineering 
The Pennsylvania State University 
University Park, Pennsylvania 


This research is motivated by the improved ballistic performance of large- 
caliber guns using stick propellant charges. A comprehensive theoretical model 
for predicting the flame-spreading, combustion , and grain deformation phenomena 
of long unslotted stick propellants is presented. The formulation is based upon a 
combined Euler lan-Lagrangian approach to simulate special characteristics of the 
two-phase con±>ustion processes in a cartridge loaded with a bundle of sticks. The 
model considers five separate regions consisting of the internal perforation, the 
solid phase , the external interstitial gas phase, and two lumped parameter re- 
gions at either end of the stick bundle. For the external gas— phase region, a set 
of transient one-dimensional fluid-dynamic equations using the Eulerian approach 
is obtained; governing equations for the stick propellants are formulated using 
the Lagrangian approach. The motion of a representative stick is derived by con- 
sidering the forces acting on the entire propellant stick. The instantaneous 
temperature and stress fields in the stick propellant are modeled by considering 
the transient ax i symmetric heat -conduct ion equation and dynamic structural analy- 
sis. For the internal perforation region, a set of one-dimensional transient 
fluid-dynamic equations is formulated with a coordinate system attached to the 
moving stick. Major distinctions between the present and the conventional formu- 
lations for interior ballistic simulation are delineated. 


I. Introduction 

Recently, there has been an increasing inter- 
est in the use of stick propellant charges in 
large-caliber gun systems. Stick propellants offer 
many advantages over conventional randomly packed 
multi-perforated granular propellant charges. The 
regular geometry of stick propellants allows a 
higher loading density, flexibility in charge 
design, and easier charge loading. The higher 
charge density is preferable for low vulnerability 
ammunition (LOVA) propellants, which require a 
higher propellant mass to produce an equivalent 
performance. It has been observed by Robbins et 
al. (1,2) that flow resistance through the charge 
of a stick propellant bundle is lower than that 
through packed beds of granular propellants, thus 
enabling faster and more reproducible flame spread- 
ing through the stick propellant charge. The lower 
flow resistance also reduces considerably the phe- 
nomena of high pressure gradients and severe 
pressure waves in the gun tube, as indicated by 
Minor (3) . 

A number of studies (1-10) on various aspects 
of stick propellant combustion have been reported 
to date. The NOVA code, developed by Gough (11,12) 
for ballistic performance of granular charges, was 
used with some modifications (1,4) to predict the 
performance of a stick propellant charge. Results 
obtained were in good agreement with experimental 
data for multiperf orated granular NACO propellants 
and single-perforated slotted stick bundles (1) . 
However, the same is not true for the case of 
single-perforated, unslotted stick propellants. 

The structure mechanics consideration in the 
continuum modeling of unslotted stick charge com- 
bustion of the modified NOVA code (4) is rather 
crude, due to the application of a steady-state re- 
lationship between radial and hoop stresses and 
internal and external pressures. Although pressure 


distributions in the internal perforation were cal- 
culated, only the external pressure was used in 
evaluating the axial stress component, which is in 
turn related to the intra-granular stress. Grain 
deformation and fracture of unslotted long sticks 
are mainly due to the radial expansion and attain- 
ment of a critical hoop stress. 

Rupture of stick propellants during combustion 
in a gun barrel was observed by Robbins and Horst 
(2) . Grain fracture can lead to high peak pressu- 
res due to increase in total burning surface area. 
The pressure difference across the web of a stick 
propellant can cause the grain to deform prior to 
fracture and alter the flow-channel width and the 
distance between opposite burning surfaces, thereby 
changing the combustion process. Hence, the pred- 
iction of grain deformation and rupture should be 
given due importance in the overall interior bal- 
listic cycle. 

The proposed model is based on a relatively 
new method, previously applied to two-phase react- 
ing flow problems such as spray combustion of 
liquid droplets (14). In the separated -flow ap- 
proach, the continuous phase (gas phase) is treated 
by a Eulerian approach, while the condensed phase 
media (stick propellants) are divided into several 
representative groups and tracked, using a Lagran- 
gian approach, as they move in the continuous 
phase. The interstitial gas-phase region in the 
present problem is treated in a similar fashion to 
the continuous phase in spray combustion, while the 
internal gas-phase regions in the grain perfora- 
tions are treated separately. The flame -spreading 
and combustion phenomena inside the perforation are 
similar to those in a cylindrical side-burning 
rocket motor grain. The ignition transient analy- 
sis developed by Peretz et al. (15) is therefore 
adapted to model the flame spreading and combustion 
processes in the perforation. This two-phase sep- 
arated flow approach is convenient because the 
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stick propellants have identical dimensions jnd 
syrrmetry scout their own axes as well as a central 
axis. The number of representative sticks in a 
stick propellant bundle is relatively low, there- 
fore making it possible to study some of the 
detailed flow and combustion phenomena of these 
representative stick propellant grains. 

All of the gun interior ballistic codes and 
analyses developed so far do not consider the dis- 
tribution of chemical species (16-18). The heat 
release is assumed to occur at the same location as 
that of pyrolysis of the propellant. However, for 
LOVA or other types of modern propellants, seme 
chemical species could be pyrolyzed from the pro- 
pellant surface, be carried by the flow during the 
ignition transient interval, and then react at some 
downstream location. Therefore, in general, the 
heat release could occur at a different position 
than the site of initial pyrolysis. These complic- 
ated phenomena, which were simplified in previous 
analyses, are included in the present model. 

Objectives of the current research are to for- 
mulate a combined Euler lan-Lagrangi an model for the 
combustion of a bundle of stick propellants inside 
a gun chamber. It is intended to cover many as- 
pects of realistic simulation of stick propellant 
combustion characteristics so that the model is ca- 
pable of predicting the phenomena of ignition, 
flame spreading, and combustion of stick propel- 
lants. Development of the model is also intended 
to help explain the experimental observations from 
test firings being conducted by the authors in a 
simulated gun system. The data to be obtained will 
be used for model validation. 

II. Method of Approach 

1. Physical Model 

Figure 1 shows a bundle of stick propellants 
loaded in the combustion chamber of a gun. For the 
present model, only unslotted long stick propellant 
grains are considered. In the theoretical formul- 
ation, the combustion chamber is divided into five 
separate regions: 1) lumped parameter region near 
the base pad; 2) internal perforation region; 3) 
external interstitial gas-phase region; 4) solid 
propellant region; and 5) lumped parameter region 
near the base of the projectile. Although the 
solid propellant region consists of many stick pro- 
pellants, only a few representative sticks are 
required to be modeled due to the similarity of the 
sticks in the same family. Each region has a sep- 
arate set of governing equations which are coupled 
through various boundary conditions. The erosive 
burning of the stick propellants under cross-flow 
conditions is also taken into account. Since frac- 
ture phenomena of stick propellants under dynamic 
loading conditions are under investigation, the 
present model is limited to the time perird before 
one or more stick propellants rupture as a result 
of high pressure differential across the web of the 
stick propellant. After the onset of rupture, a 
stick propellant could become partially slotted, 
broken, and/or highly deformed. It is then diffi- 
cult to distinguish the internal versus external 
surface. These phenomena are beyond the scope of 
the present model. 

2. Mathematical Formulation 
A. Basic Assumptions: 
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Fig.] Stick propellant Charge in a Cartridge 
of a Large Caliber Gun. 


(1) All the stick propellants have the same 
initial geometry and physical conditions. Further- 
more, it is possible to divide the stick bundle 
into a few families so that the calculations can be 
performed for only a few representative sticks. To 
simplify the mathematical formulation, the combus- 
tion of sticks in a bundle is represented by a 
single stick. The mathematical format can be read- 
ily extended to multiple families of sticks. 

(2) Assumptions used in the gas-phase regions 

are: 

a. no body forces; 

b. bulk viscosity ji r is negligible; 

c. Soret and CMfour effects are negligi- 
ble ; 

d. gases obey Noble -Abel equation of 
state; 

e. all binary diffusion coefficients are 
equal; 

f. Fick's law of diffusion is valid; 

g. flow is one-dimensional transient 
(properties are uniform in r and 9 di- 
rections) ; 

h. total flow area in the external gas- 
phase region is assumed to be M 
[Dupuit Forchheimer hypothesis (19) f; 
and 

i. turbulence correlations of flow proper- 
ties in the axial direction are consid- 
ered negligible in comparison with the 
product of their mean flow properties. 
However, turbulence effects in the 
transverse direction are embedded in 
the empirical correlations. 

(3) Assumptions used for the solid propellant 

are : 

a. burning in the inner and outer surfaces 
of a stick is axisynmetric; 

b. density of the solid is constant; 

c. any subsurface heat release occurs very 
close to the surface and therefore can 
be lumped onto the surface; 

d. torsion and rotation are negligible; 

e. each stick propellant is locally axi- 
synmetric; 

f. burning at the end surfaces is uniform; 

g. end surfaces are perpendicular to the 
axis of the stick; and 

h. propellant material behaves as a linear 
viscoelastic material in shear and 
elastic material in bulk deformation. 

B. 'Overall Structure of the Mathematical 
Model 
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In order >:o give an overview of the structure 
tne mathematical model, Pig. 2 was constructed. 
It shows the major components and subcomponents of 
the formulation for each region in the physical 
model. To simplify the diagram, the cross -coup ling 
‘mes between each region are not shown. However, 
it is important to note that the physicochemical 
processes in these regions are closely coupled. 
Mathematical representations of each component and 
subcomponent are given in the following sections. 

C. Governing Equations for the External 
Gas-Phase Region 

The external gas-phase region occupies the 
whole interstitial space in the cartridge, exclud- 
ing the shaded region shown in Fig. 3. By using 
the control volune analysis (reference to Eulenan 
coordinate) , the following governing equations are 
obtained . 
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Fig. 3 Schematic Diagram Showing the Mass 
Fluxes Entering and Leaving the Exter- 
nal Gas-Phase Region (Unshaded and 
Multiply Connected Region) . 


C >ncinuity Equation: 
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where 4> e is the void fraction of the external gas- 
phase region and is defined as 


1 - 


r 2 A °j 


( 2 ) 


Momentum Equation: 
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where 7 sin 3 0 is the horizontal component of the 
gastif ication velocity normal to the solid propel- 
lant surface. The second term in the right-hand 
side represents the momentum transfer from solid 
propellant to the external gas-phase region due to 
burning . 


Energy Equation: 
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vtere h. is the enthalpy of the ith species defined 
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( 5) 


^ the sum of the internal and kinetic 
energy of the gas phase in the control volume. The 
energy transfer due to molecular species diffusion 
has been neglected because of a high degree of tur- 
du lence « 
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where the source 
contributions cause 
phase reactions, i.e. 


term (tu^) consists of 
by surface pyrolysis and gas- 
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where Y^ represents the mass fraction of the ith 
species pyrolyzed from the solid propellant before 
mixing with ambient gases. Following the flame 
model proposed by WU et al. (20) in their study of 
erosive burning of homogeneous propellants, the 
solid propellant pyrolyzes into three groups of 
species: 


Solid Propellant Oxidizer-rich gases (0) 

+■ 

Fuel-rich gases (F) 

+ 

First group of species 
with delayed reaction ; 

( 9) 

Utider low cross flow conditions, the flame struc- 
ture adjacent to a burning homogeneous propellant 
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surface exposed to a Large cavity can be assumed to 
nave three stages, as shown in Fig. 4a. This flame 
structure is based upon the erosive ourmng study 
in a rocket motor by Wu et al. (20) . In the case 
of Large-caliber guns densely Loaded with stick 
propellants , the void spaces adjacent to burning 
surfaces are relatively small/ -and the species py- 
rolyzed from the surface can be entrained by the 
high-veiocity gases flowing along the axis. 'The 
heat release in the final flame generated by the 
chemical reaction of pyrolyzed species from a spe- 
cific location occurs at a downstream location, as 
shown in Fig. 4b. To determine the gas phase reac- 
tion rate and heat-release rate, the same chemical 
reaction mechanism, proposed in Ref. 20, is adop- 
ted. This mechanism can be represented by three 
overall chemical steps: 
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Oxidizer-rich gases (0) can be regarded as 
species, fuel-rich gases (F) as a group of aldehy- 
des (CH~0) f and other species generated by surface 
oyrolysis as one group of delayed reaction species 
(DRl) . Based upon the work of Fifer (21) and 
Kubota (22), chemical reaction in the Fizz zone is 
Largely due to the reaction involving reduction of 
NO to NO. After the delay in the dark zone, the 
reactions in the final flame can be assumed to 
occur at- high activation energies, associated with 
reactions (11) and (12) to form the final products 
from DRl and the second group of delayed reaction 
species (DR2) . In the dark and final flame zones, 
chemical reactions result in oxidation of CO, and 
perhaps some H 2 by NO. More detailed discussions 
of kinetic parameters and mechanisms are given in 
Ref. 20-22. 


It is useful to point out that the chemical 
reaction mechanism proposed in Ref. 20 can be fol- 
lowed for highly turbulent cross-flow situations. 
The 0 and F species pyrolyzed from the solid pro- 
pellant originate at the same place, and flow 
together in a torturous path to form DR2 and final 
products (P) . The delayed reaction species of DRl 
can also be considered to flow together in the 
process to form product species P. Gases in the 
cross-flow can transfer heat to these species and 
alter the delay times required to form product spe- 
cies P. The stoichiometric coefficients v p , v n , 
' , and ^ _ are therefore functions of prope r - 
lanc ingredients only. The method for calculating 
these parameters is based upon the original molecu- 
lar structure of the propellant and is discussed in 
Ref. 20. 


In Ref. 20, the rate of production of species 
"i" was based upon the chemical reaction rate as 
well as the eddy-break-up rate (23) controlled by 
the turbulence intensity and concentration gra- 
dient. Turbulence intensity is so high in the gun 
situation because of base-pad ignition and combus- 
tion of propellants that the species diffusion term 
can be regarded as extremely short* Therefore, the 
rate of consumption or production of species is de- 
termined solely from chemical reaction time. 


The rate of chemical reaction of species F, 
DRl, 0, and DR2 are given below in the same form as 
m Ref. 20. 
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Fig. 4 Flame Structure of a Homogeneous Solid 
Propellant (a) Two-dimensional Struc- 
ture under Low Cross-Flow 

Velocities, (b) Distorted Flame under 
Extermely High Cross-Flow Velocities 
(Nearly One-Dimensional Structure) . 


Equation of State: 

The Noble-Abel dense gas law is used. 

P (_i b) = RT (17 > 

e S e 

The initial and boundary conditions as well as em- 
pirical correlations for the external gas-phase 
region are given in a later section. 

D. Governing Equations for the Stick 
Propellant 

The following governing equations for a repre- 
sentative stick are derived based upon the 
Lagrangian coordinate. 

1) Mass and Momentum Equations: 

The instantaneous mass of the representative 
stick can be calculated from 
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L( t) / 2 


2) Transient Heat Conduction Equation: 
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The instantaneous values of stick propellant 
length, L(t), and the local inner and outer radii, 
r. (t,£) and r (t,£), ace determined by integrating 
tAe following°first order differential equations. 
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To determine the instantaneous temperature 
distribution in the stick propellant, a heat con- 
duction equation must be considered. The equation, 
which takes into account the subsurface radiation 
absorption for translucent propellants, has the 
following form. 
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where V and V represent the radial velocities 
of the sr lnner aRS°outer surfaces of the stick pro- 
pellant due to mechanical deformation with respect 
to the centerline of the stick. V gRB and re- 

present the rate of mechanical deformation or the 
right and left boundary surfaces with respect to 
the geometric center of the stick. The equation of 
motion is formulated according to the following 
momentum balance principle. 

d (M u ) 

— — — «= £ + (net rate of momentum flux flowing 
into the control volume encompass- 
ing the stick propellant) 

( 22 ) 

which gives (see Fig. 5j 


where E, i| the black-body arussive power evaluated 
as oT , and "a" represents the flux absorption 
coefficient of the propellant. This equation is 
based i£>on a two-flux model which assumes that the 
radiation fluxes are dominant in the radial direc- 
tions (inward and outward). The source terms 
represent the net rate of energy absorbed due to 
radiant energy fluxes. The outward and inward 

radiant fluxes, I and J , can be determined from 
r r 

the following flux -transport equations. 

d(rl ) 1 

— = -(s+a) rl + arE + J + J sr(I +J ) (25) 

dr r d r c 

d(r J ) i 

— - — — = (s+a)rJ - arE + J - y sr(I +J ) (26) 

dr r d r ^ l t 

These two equations were used by Gosman and Lock- 
wood . (24) 
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Where F pw re P cesents the net force acting between 
the stick propellant and chamber wall, and F is 
the force between adjacent propellants. pp 


3) Dynamic Structural Analysis 

As a result of the different ignition and 
flame spreading processes occuring in the internal 
perforation and exterral gas-phase regions, a pres- 
sure differential exists across the web of a 
representative stick. A finite element analysis is 
needed to compute the resulting transient viscoe- 
lastic deformation of the stick propellant and to 
predict the attainment of a critical condition for 
grain fracture. Regression of the boundary as a 
result of pyrolysis and burning should also be 
taken into account. Furthermore, the mechanical 
properties of the stick propellant must be specif- 
ied. 


The propellant material can be considered as 
linear viscoelastic in shear and elastic in bulk 
deformation. This is a conmonly accepted practice 
for solid propellant (25, 26). The elastic bulk 
behavior is assumed to follow 
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Fig. 5 Momentum fluxes and Pressure Forces 
Acting on the Entire Surface of a 
Single Perforated Stick Propellant. 


where K is the bulk modulus. The deviator ic beha- 
vior is taken as 

( C 3e (f) 

S ij ” G i (t - t ’> -dt' ( 28 ) 

J O 

where the shear relaxation modulus G. (t) is assumed 
to be of the form 

G ( t) « G + (G -G )e" gt (29) 
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where is the long-time shear modulus , G q is the 
short-time shear modulus, and S is the decay con- 
stant. Since a closed-form solution for the 
axisynmetric dynamic problem is not possible, a 
well-established finite element code "HONDO- I I" 
(26) is employed for computations of grain deform- 
ation. 1 he code utilizes the principle of virtual 
work for the solution. It states that at all the 
points along the path of motion, the differential 
virtual work, 6 it r must vanish for all variations 
satisfying the imposed displacement boundary condi- 
tions (26). Stt is defined as 
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where t* is surface traction and a is Cauchy 
stress tensor. The stick propellant is divided 
into a number of elements over the cross-sectional 
area of the web. The HONDO- I I code uses four mode 
bilinear isoparametric elements. The basic equ- 
ation of motion, viz., the minimization of virtual 
work, is then considered for each of the nodes in 
the stick propellant. Thus, the equation of motion 
for a node becomes 
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and <t> are the bilinear interpolation functions. 
In Eq. (31), W is the total number of elements sur- 
rounding the node in question. The calculations 
are carried out on an element -by -element process to 
get the final equations of motion. The time in- 
tegration of these equations gives the positions of 
nodes at the new time step. A central difference 
method is used for time integration in the HONDO-II 
code. 


E. Governing Equations for the Internal 
Gas-Phase Region 

Conservation equations for the control volume 
in the perforation of a stick propellant are simi- 
lar to those given in Ref. 15. 
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Momentum Conservation Equation: 
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Energy Conservation Equation: 


The energy conservation equation written for 
the stored total energy (internal and kinetic) per 

unit mass, E . , is 
gi' 
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Species Continuity Equation: 


The species continuity equation for the inter- 
nal gas-phase is similar to that for the external 
gas-phase, with the void fraction taken as one and 
the absolute velocity replaced by the relative ve- 
locity. 
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where the source term (u. )., similar to ( w . ) in 
Eq. (7) , consists of contributions due to surface 
pyrolysis and gas-phase reactions. The flame model 
for the internal gas-phase region is the same as 
that given earlier for the external gas-phase 
region. 


F. Heat Losses to the Walls of the Combustion 
Chamber 


In order to consider heat losses from the com- 
bustion zone to the cartridge chamber, the gun 
tube, and the projectile, temperature profiles in 
these metal components are required. For a test rig 
with a blowout diphragm and short barrel, the tran- 
sient ignition and combustion phenomena occur in an 
extremely short time interval and the heat loss to 
the surrounding walls can be considered negligible. 
However, if the simulation is made for the full 
ballistic cycle occurring in a long gun barrel, the 
loss to the walls must be considered. The heat con- 
duction to the gun tube can be given as 
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where £ is the distance measured from the center of 
the stick, and V • is the velocity of the gas rela- 
tive to the velocity of the stick propellant. V • 
is related to the absolute velocity u ^ by 


The transient heat conduction equation for the 
wall on the breech end will follow the one- 
dimensional form: 
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2). Boundary Conditions for the External and 
(38) Internal Gas-Phase Regions 


The equation for the projectile is: 


D(C T ) . 3t 

w w d ( w 

P w DC 3 3z 3z 


(39) 


The solution of these equations are coupled to 
the external gas-phase region through the boundary 
conditions on the wall surfaces. 


For simplicity in expressing the boundary 
treatments, only one representative group of stick 
propellants is considered, then each stick propel- 
lant has the identical instantaneous velocity and 
length. Therefore, the governing equations for the 
left lumped -parameter region (see Fig. 6) or called 
left control volume (LCV) can be readily derived 
and written as follows. 


G. Initial and Boundary Conditions 

For each of the above governig equations, 
there is a set of boundary and/or initial condi- 
tions required to complete the formulation. The 
initial conditions can be specified readily, since 
the gas and stick propellant velocities are zero 
and the pressure and temperature of gas in the car- 
tridge is at room conditions. The stick propellant 
charge with known geometry and its surrounding 
chamber are also at room temperature. Since the 
amount of .air in the initial loading of the car- 
tridge is extremely small in comparison with the 
gases generated from combustion, it can be treated 
as any one of the five species discussed above, in 
view of the fact that air contains oxygen, it is 
treated as oxidizer-rich gas (0) . 


Mass balance: 


LCV 


dt 


= - P, 


(u +r ) 
3 b LB 


8 lcv z l 

r 


l-V ~ 


6 LCV 

<T) 

P . LCV 
ign j 


m t -m -m 
ign oute out! 


UR c \CV Z L 


(44) 


where is the length of the left lumped -parameter 
region and can be calculated from 


Z L = Z L + 
o 


and 


(u +r )dt + 
3 b LB 




n 2 
™ P. 
o c ign 


dt (45) 


A number of important boundary conditions are 
given in the following. 

1). Boundary Conditions for the stick Pro- 
pellant 


The boundary condition for Eq . (24 ) at outer 
surface of the stick propellant pan be written as 


3T 


s9 r 


r h ( r -T ) + p r (Q ) + 

r 0 c g e se s be s chem ^rad'e 




+ P r, (C -C )T + (Q ) (40) 

s be s pge se V CPIP e 


where ((^p.p) re P resen ts the rate of energy input 
due the deposit of condensed ph$se igniter products 
onto the external surface and (q ) the net ra- 
diative heat flux to the surface can be expressed 
as 


7TR p [U - (U +r ) ; 


oute e c pg ge s b 


LB 


i = N7Tr i p g i [,, gi- ( V r b 


out i 


Momentum balance: 


LB 


Z - -L/2 
Z « -L/2 


(46) 

(47) 
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Similarly, the boundary condition for the inner 
surface of the stick propellant can be written as 


- k 


9T 

s 

s 3r 


, = h (T -T ) + p r. . (Q ) + (q ,) 

+ c gi si s hi H s chem rad 

r i 1 


+ P s r bi (C s - C pgi )T si + <Wi (42) 


where 

<Oi “ 


+ £ J 
s r 


t - e E, T 
+ s b si 


(43) 


The pressure distributions along the internal per- 
foration and external interstitial regions solved 
from the gas-phase equations are used as boundary 
conditions ■ for the solid-phase dynamic structral 
analysis. 


<P„ 
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Fig. 6 Control Volume Considered for Deriving 
Governing Equations in the Left 
Lumped-Parameter Region. 
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Energy balance: 


region remains the same as that in the first case. 


dT 


LCV 


dt (C vV L CV 
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o 
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> (49) 


Equation of state: 


RT 


& LCV 


LCV 


- b 


(50) 


The void fraction in the left limped -parameter 
region, can be evaluated from the continuity 
equation or the condensed -phase products of the ig- 
niter which is given as follows. 


d<t>. 


LCV 


dt 


m . - ( 1 - 


(l-b ) f u 4r 4 A&0 — ] . 

V acv ; 1 s b n 2 j 
LB TTR p, 

c ign 2 

) [p A$ u +p Ntit . u 


len ' S ' LCV lgn s e LB 8 e LB 1 8l LB 

11 R 2 p. 

C 


(51) 


The boundary values of velocity, density, 
pressure, and temperature on the left hand side of 
both internal and external gas-phase regions can be 
solved from a number of relationships coupled with 
Eqs. (44) and ( 48 ) — (51 ) . Depending upon the flow di- 
rections at left boundary of the internal and 
external regions, three different cases are identi- 
fied. In the first case, the flow velocities at 
stick perforation and the interstitial regions are 
both in the positive z direction, it can be assumed 
that the pressure at LCV is equal to those at the 
left boundary of the stick perforation and the ex- 
ternal gas-phase region. Also, the stagnation 
temperature of the gas in LCV can be assumed to be 
the same as those on the left boundary. Oie com- 
patibility relation along the characteristic line 
in each gas-phase region is used together with the 
aforementioned equations for solving the flow prop- 
erties in the LCV and the left boundary. 


In the third case, the flow velocities at the 
left boundray of both the perforation region and 
the interstitial regions are negative. Therefore, 
two additional compatibility relationships are 
available. In this case, both the internal and ex- 
ternal gas-phase boundaries are treated in a 
similar fashion as that for the internal perfor- 
ation region in the second case. 

The void region between the right end of the 
stick propellants and the projectile base is a 
lumped -parameter region, which could contain either 
a base pad with powders of igniter material or an 
ullage space. In the case of base pad, a set of 
equations similar to Eqs. (44)- (51) are derived. The 
right control volume (RCV) is considered to move 
with the projectile and the left surface of RCV 
coincides with the right boundary of the stick 
bundle. A detailed discussion on the compatibility 
relation mentioned above is given in Ref. 27-28. 

H. Empirical Correlations and Constants 

Several empirical correlations are used in the 
model. For the internal gas-phase region, these 
comprise of the correlations for the burning rate 
law including erosive burning effect, the convec- 
tive heat-transfer coefficients, and viscous drag 
coefficient. The latter two correlations are same 
as those used in Ref. 15. New erosive burning cor- 
relation being developed by authors and coworkers 
in parallel to this study will be used. For the ex- 
ternal gas-phase region, correlations for the 
convective heat-transfer coefficient, and viscous 
drag coefficient between the gas-phase and solid 
surfaces are needed. These correlations are similar 
to those for the internal perforation region except 
the fact that they are based on a hydraulic diame- 
ter D defined as 

H 

2 2 
2 ( R - Nr ) 

rx - C O 


c o 

The burning rate law for the external region 
will be the same as that for the internal perfor- 
ation region. 

For the stick propellants, correlations for 
F and F , the ignition criterion, and mechanical 
aB8 therm&f properties are needed. At the present 
time, not all of this information is available; es- 
pecially F and F need to be characterized. 

PP pw 

4. Summary of Differences Between the 
Present and Conventional Formulation 

In order to bring out the major differences 
between the present formulation and the conventio- 
nal interior ballistic predictive models (4), a 
summary table is given below. 


In the second case, the flow velocity at the 
left opening if the stick perforation is negative, 
but the bulk velocity at left boundary of the in- 
terstitial void region is still positive. In this 
case, one additional compatibility relationship can 
be used .for the stick perforation region. Only the 
pressure at the opening of the stick perforation is 
considered to be equal to that in the LCV. The 
treatment for the left boundary of the interstitial 


As sunmerized by Robbins and Einstein (29), 
there are differences between the measured and cal- 
culated pressures from the NOVA code or its 
extensions for long unslotted stick propellants. 
Also, the measured muzzle velocities are higher 
than those calculated for slotted stick propel- 
lants. A number of improvements and considerations 
suggested in the workshop (29) are incorporated in 
the present model. 


Of 
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Differences Between the Present Formulation ana the 
Conventional Interior Ballistic Formulation 





Subject under 
Consideration 

Present Formulation 

Conventiai formulation (4) 

Typical grain 
configuration 

* Simulation of a number of 
typical full-length grains 
in a bundle of stick prop- 
ellants. 

* Each stick is modeled as 
separate tube with deform- 
able and combustible walls. 

* Simulation of an average 
grain in a spatial location 
along a packed bed of gran- 
ular propellants, 

* Each bundle is modeled as a 
continuum characterized by 
the velocity and stress in 
the sticks. 

Averaging of 
flow 

properties 

* The external flow propert- 
ies are averaged over the 
cross-sectional flow area 
of interstitial voids, 
while the Internal flow 
properties are averaged 
over the flow area of each 
stick propellant. 

* Flow properties in both ext- 
ernal and internal regions 
are averaged over their res- 
pective flow areas. 

1 

Grain 

deformation 
and fracture 

* Simulated by the unbalanced 
pressure forces between the 
internal perforation and 
external interstitial void 
region. 

* Linear viscoelastic consti- 
tutive law is used. 

* Employs dynamic finite- 
element structure mechanics 
computational code. 

* The process of grain defo- 
rmation and fracture are 
not addressed, except the 
longitudinal stresses are 
considered in the solid- 
phase momentum equation. 

* Linear elastic constitutive 
law is used. 

* Employs a steady-state rel- 
ationship between stresses 
(radial and hoop) and press- 
ures (internal and external) J 

Gram 

displacement 

and 

acceleration 

* The kinematics of the full- 
length grain is determined 
from the summation of all 
forces exerted on the grain. 

* The bulk properties of the j 

grains are determined from • 

local momentum balance. 

1 

Radiative heat 
transfer 

; * Subsurface radiation pene- 
tration is allowed and tre- 
ated by a two-flux model. 

: 

* No subsurface radiation j 

penetration is considered. 1 

i 

Type of 
; formulation 
and frame of 
reference 

* Kinematics and grain defor- 
mation are formulated by 
following the stick (Lagra- 
ngian approach) while the 
gas-phase properties for 
internal and external regi- 
ons are determined from a 
fixed frame of reference 
{Euler ian approach) . 

* Both the gas-phase and solid 
-phase properties are all 
determined from the conserv- 
ation equations formulated 
based upon a fixed frame of 
reference (Euler ian 
approach) . 

Species 
distribution 
and location of 

* Five groups of species are 
considered. 

* Gas phase is made of combus- 
tion products from ignition 
and propellant. 

heat release 

* Heat release does not have 
to occur at the site of 
pyrolysis. 

* Heat release occurs locally 
at the site of pyrolysis. 

[ 
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Letter Symbols 


Nomenclature 


a flux model absorption coefficient, m 
A cross-sectional area of the gun £grrel, rn 

preexponential factor of the i species, 
m /kmol-s ^ 

Ap cross-sectional area of the perforation, m 

A g specific area of the external surfaces of 
the stick propellant, m 

b covolume, of the Noble -Abel equation of 
state, m /kg 

0^ constant pressure specific heat, J/kg-K 

C s specific heat of stick propellant, JAg-K 

C v constant volute specific heat, JAg-K 

DRl group of species pyrolyzed from propellant 
surface having delayed reactions 
DR2 delayed reaction species generated from 0 
and F species 

D v viscous drag force per unit area, N/nr 

P binary diffusion coefficient, m 2 /s 

deviatonc strain 

E total stored energy (internal plus kinetic) 
per unit mass, JAg 
activation energy, J/kmol 

4 2 

black-body emissive power, = ; T , J/m -s 

f^ body-force per unit volume in k th direction, 
N/rn 

F fuel rich species pyrolyzed from propellant 

F^ external force exerted on the solid propel- 
lant in the axial direction, N 
h specific enthalpy, JAg 

h c convective heat-transfer coefficient, W/iri -K 

fu specific enthalpy of the species, JAg 

h fc total heat-transfer coefficient, W/m -K 

0 th 

Ahf i standard enthalpy of formation of the i 

' species, JAg 


1^ outward radiation flux in the positive 
1 radial direction, W/rn 

J^ inward radiation flux in the negative radial 

2 

direction, W/m 

k thermal conductivity, W/m-K 

K bulk modulus of the propellant material, 

N/nr 

Mg total mass of gas in control volute, kg 

M g instantaneous total mass of a single stick 
propellant, kg 

0 oxidizer species pyrolyzed from propellant 
P pressure, N/rn 

P final product species 

perimeter of the internal perforation, m 

» 

radiative heat fulx per unit time ^bsorbed 
by the solid propellant surface, W/nri 
Q surface heat release due to pyrolysis, JAg 
chem 

* 2 
0 rate of heat loss to the tube wall, W/m 


r^ propellant burning rate, m/s 
r^ inner radius of the perforation, m 

r Q outer radius of the stick propellant, m 


R gas constant, J/kg-K 

R radius of the combustion chamber, m 

c 


R^ universal gas constant, J/kmol-K 




S 

t 

T 


g 

i/ 

w. 


flux-model scattering coefficient, 1/m 
deviatonc stress tensor, N/m 

2 

surface, m 
time, s 

temperature, K 

adiabatic flame temperature of the solid 
propellant, K 
absolute velocity, m/s 

gas velocity relative to the solid propel- 
lant, m/s 3 

volume, m J ^ 

molecular weight of ££je i species, kg/kmol 

coordinate axis in k direction, m 

acceleration in k^ direction, m/s 2 


5x^ m deformation gradient tensor 

mass fraction of i* 1 * 1 species, l could 
represent F, 0, DRl, DR 2, or P 
z axial coordinate in Eulecian coordinate 
system, m 


Greek Symbols 
dilatory strain 

surface emissivity of solid propellant 

angle measured ccw from axis to the 

tangent to the inner surface of the stick 
propellant, rad 

0^ angle measured ccw from axis to the tan- 
gent to the outer surface of the stick 
propellant, rad 

p dynamic viscosity of the gas, N-s/m 2 
number of kmoles of l species 

£ Lagrangian axial coordinate, m 
tt virtual work, J 

p density, kg/rn 

Stefan-Boltzmann constant, W/m-K 4 
a stress tensor, N/nr 

2 

t viscous shear stress, N/m 

T k surface traction in k th direction, N/m 2 

2 

t w tube wall shear stress, N/m 

2 

t £-£■ normal viscous stress ♦ N/m 

mass fraction of combustion product in the 

condensed phase from th<= igniter 

rate of production of l species, kg/m -s 

Subscripts 

e external interstitial region 

g gas-phase, internal or external region 

i internal perforation region 
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- : n 

igniter 




left boundary of 

the stick propel 

.ant bundle 

LCV 

left lumped -parameter region in 

the car- 


triage 



RB 

ngnt boundary 
bundle 

of the stick 

propellant 


solid propellant 

or surface 
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Editor's Narrative: 


Round-Table Discussion: Propulsion Applications of Mixing and 
Demixing Processes of Multiphase Flows 

Several points were made either explicitly or implicitly in 
the presentations and discussions on the second day concerning 
mixing processes in multiphase systems and the modeling of these 
systems. First, there are situations where transport in turbu- 
lent flows can be in a direction other than along the gradient of 
the property being transported. This comes about through cou- 
pling of density inhomogenieties with pressure gradients, the 
former being a consequence of the large temperature differences 
in combustion environments. Hence, in turbulent flames, some 
transport processes are "counter gradient" or more correctly, 
"non-gradient" in nature. Errors due to neglect of these pro- 
cesses are not understood, but could possibly be significant in 
some combustion systems. Modelers are in the process of learning 
to incorporate these ideas into practical codes. 

Second, the weakest element in modeling of turbulent multi- 
phase flows is probably the turbulence modeling itself. The 
treatment of the coupling between the turbulence and the non- 
continuously distributed dispersed phase is also a critical area 
and includes the question of whether, for numerical purposes, to 
treat the dispersed particulate (or droplet) phase as a continuum 
or as an assembly of separate particles. 


PRECEDIN': 


1 h * 1 V 
Zj 
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Third, it is possible to treat a very complex problem in 
great detail by using appropriate levels of approximation and by 
coupling separate regions together which are best treated with 
different models or with different levels of approximation. 

Fourth, while continuum approaches to multiphase flow calcu- 
lations are not particularly illuminating in revealing details of 
the physical processes at work (due to the necessity of averaging 
quantities before solving the governing equations) , there are 
cases where they can provide useful results for a given problem. 
The two-fluid approach works best where the dispersed particulate 
phase is monodisperse , and where the flows are non— reacting . 

Sixth, advantages and disadvantages of Favre averaging 
approaches were discussed. An advantage is that the conservation 
equations for variable density with Favre averaging are much like 
the standard constant density equations. Disadvantages are that 
there is some difficulty in obtaining molecular terms and com- 
puted and measured quantities are more difficult to compare. 

Finally, physical situations for which multiphase flow mod- 
els are appropriate generally contain a very broad range of 
length and time scales within the same problem. This means that 
when computational approaches to these problems are used, some 
levels of approximation will need to be made in order to deal 
with these large ranges of scale. The point in the analysis at 
which approximation is made is a key difference between a two- 
fluid model and an Eulerian— Lagrangian model of a dispersed par- 
ticulate phase in a fluid phase. 
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